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(57) Abstract 

A technique is described for determining the ef- 
fectiveness of medical therapy and dosage formulations 
by analysing dot spectrograms representative of quan- 
tized hybridization activity in biological samples, such 
as DNA, RNA or other protein biomolecular array sam- 
ples, taken at different sampling times from a patient 
undergoing the treatment. The technique directly lends 
itself to disease progression analysis based on markers 
such as viral load. In accordance with the technique, 
a viral diffusion curve associated with a therapy of in- 
terest is generated and each dot spectrogram is then 
mapped to the viral diffusion curve using fractal fil- 
tering to yield a filtered viral diffusion curve for each 
sample. A degree of convergence between the filtered 
viral diffusion curves is determined. Then, a determi- 
nation is made as to whether the therapy of interest 
has been effective by determining whether the degree 
of convergence increases from one sample to another, 
with an increase in the degree of convergence being 
representative of a lack of effectiveness of the therapy 
of interest. 
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5 METHOD AND APPARATUS FOR EXPONENTIALLY CONVERGENT THERAPY 

EFFECTIVENES MONITORING USING DNA MICROARRAY BASED VIRAL LOAD 

MEASUREMENTS 
Field of the Invention 

The invention generally relates to techniques for monitoring the effectiveness of medical 

10 therapies and dosage formulations, and in particular to techniques for monitoring therapy 
effectiveness using viral load measurements. 

Background of the Invention 
It is often desirable to determine the effectiveness of therapies, such as those directed against 
viral infections, including therapies involving individual drugs, combinations of drugs, or other 

1 5 related therapies. One conventional technique for monitoring the effectiveness of a viral infection 
therapy is to measure and track a viral load associated with the viralinfection, wherein the viral 
load is a measurement of a number of copies of the virus within a given quantity of blood, such per 
milliliter of blood. The therapy is deemed effective if the viral load is decreased as a result of the 
therapy. A determination of whether any particular therapy is effective is helpful in determining the 

20 appropriate therapy for a particular patient and also for determining whether a particular therapy is 
effective for an entire class of patients. The latter is typically necessary in order to obtain FDA 
approval of any new drug or medical device therapy. Viral load monitoring is also useful for 
research purposes such as for assessing the effectiveness of new antiviral compounds determine, for 
example, whether it is useful to continue developing particular antiviral compounds or to attempt to 

25 gain FDA market approval. 

A test to determine the viral load can be done with blood drawn from T-cells or from other 
standard sources. The viral load is typically reported either as an absolute number, i.e., the number 
of virus particles per milliliter of blood or on a logarithmic scale. Likewise, decreases in viral load 
are reported in absolute numbers, logarithmic scales, or as percentages. 

30 It should be noted though that a viral load captures only a fraction of the total virus in the 

body of the patient, i.e., it tracks only the quantity of circulating vims. However, viral load is an 
important clinical marker because the quantity of circulating virus is the most important factor in 
determining disease outcome, as changes in the viral load occur prior to changes in other detectable 
factors, such as CD4 levels. Indeed, a measurement of the viral load is rapidly becoming the 

35 acceptable method for predicting clinical progression of certain diseases such as HIV. 

Insofar as HIV is concerned, HlV-progression studies have indicated a significant 
correlation between the risk of acquiring AIDS and an initial HIV baseline viral load level. In 
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5 addition to predicting the risk of disease progression, viral load testing is useful in predicting the 
risk of transmission. In this regard, infected individuals with higher viral load are more likely to 
transmit the virus than others. 

Currently, there are several different systems for monitoring viral load including quantiti ve 
polymerase chain reaction (PCR) and nucleic acid hybridization. Herein, the term viral load refers 

1 0 to any virological measurement using RNA, DNA, or p24 antigen in plasma. Note that viral RNA 
is a more sensitive marker than p24 antigen. p24 antigen has been shown to be detectable in less 
than 50% asymptomatic individuals. Moreover, levels of viral RNA rise and fall more rapidly than 
levels of CD4+ lymphocytes. Hence, changes in infection can be detected more quickly using viral 
load studies based upon viral RNA than using CD4 studies. Moreover, viral load values have to 

1 5 date proven to be an earlier and better predictor of long term patient outcome than CD4-cell counts. 
Thus, viral load determinations are rapidly becoming an important decision aid for anti-retro viral 
therapy and disease management. Viral load studies, however, have not yet completely replaced 
CD+ analysis in part because viral load only monitors the progress of the virus during infection 
whereas CD4+ analysis monitors the immune system directly. Nevertheless, even where CD4+ 

20 analysis is effective, viral load measurements can supplement information provided by the CD4 
counts. For example, an individual undergoing long term treatment may appear stable based upon 
the observation of clinical parameters and CD4 counts. However, the viral load of the patient may 
nevertheless be increasing. Hence, a measurement of the viral load can potentially assist a 
physician in determining whether to change therapy despite the appearance of long term stability 

25 based upon CD4 counts. 

Thus, viral load measurements are very useful. However, there remains considerable room 
for improvement. One problem with current viral load measurements is that the threshold level for 
detection, i.e., the nadir of detection, is about 400-500 copies per milliliter. Hence, currently, if the 
viral load is below 400-500 copies per milliliter, the virus is undetectable. The virus may 

30 nevertheless remain within the body. Indeed, considerable quantities of the virus may remain 
within the lymph system. Accordingly, it would be desirable to provide an improved method for 
measuring viral load which permits viral load levels of less than 400-500 copies per milliliter to be 
reliably detected. 

Another problem with current viral load measurement techniques is that the techniques are 
35 typically only effective for detecting exponential changes in viral loads. In other words, current 
techniques will only reliably detect circumstances wherein the viral load increases or decreases by 
an order of magnitude, such by a factor of 1 0. In other cases, viral load measurements only detect a 
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5 difference between undetectable levels of the virus and detectable levels of the virus. As can be 
appreciated, it would be highly desirable to provide an improved method for tracking changes in 
viral load which does not require an exponential change in the viral load for detection or which does 
not require a change from an undetectable level to a detectable level. Indeed, with current 
techniques, an exponential or sub-exponential change in the viral load results only in a linear 

10 change in the parameters used to measure the viral load. It would instead be highly desirable to 
provide a method for monitoring the viral load which converts a linear change in the viral load into 
an exponential change within the parameters being measured to thereby permit very slight 
variations in viral load to be reliably detected. In other words, current viral load detection 
techniques are useful only as a qualitative estimator, rather than as a quantative estimator. 

15 One reason that current viral load measurements do not reliably track small scale 

fluctuations in the actual number of viruses is that a significant uncertainty in the measurements 
often occurs. As a result, individual viral load measurements have little statistical significance and 
a relatively large number of measurements must be made before any statistically significant 
conclusions can be drawn. As can be appreciated it would be desirable to provide a viral load 

20 detection technique which can reliably measure the viral load such that the statistical error 
associated with a single viral load measurement is relatively low to permit individual viral load 
measurements to be more effectively exploited. 

Moreover, because individual viral load measurements are not particularly significant when 
using current methods, treatment decisions for individual patients based upon the viral load 

25 measurements must be based only upon long term changes or trends in the viral load resulting in a 
delay in any decision to change therapy. It would be highly desirable to provide an improved 
method for measuring and tracking viral load such that treatment decisions can be made much more 
quickly based upon short term trends of measured viral loads. 

As noted above, the current nadir of viral load detectability is at 400-500 copies of the virus 

30 per milliliter. Anything below that level is deemed to be undetectable. Currently the most 
successful and potent multi-drug therapies are able to suppress viral load below that level of 
detection in about 80-90 percent of patients. Thereafter, viral load is no longer an effective 
indicator of therapy. By providing a viral load monitoring technique which reduces the nadir of 
detectability significantly, the relative effectiveness of different multi-drug therapies can be more 

35 effectively compared. Indeed, new FDA guidelines for providing accelerated approval of a new 
drug containing regimen requires that the regimen suppress the viral load below the current nadir of 
detection in about 80 to 85 percent of cases. If the new regimen suppresses the viral load to 
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undetectable levels in less than 80 to 85 percent of the cases, the new drug will gain accelerated 
approval only if it has other redeeming qualities such as a preferable dosing regimen (such as only 
once or twice per day), a favorable side effect profile, or a favorable resistance or cross-resistance 
profile. Thus, the ability of a regimen to suppress the viral load below the level of detection is an 
important factor in FDA approval. However, because the level of detectability remains relatively 
high, full approval is currently not granted by the FDA solely based upon the ability of the regimen 
to suppress the viral load below the minimum level of detection. Rather, for full approval, the FDA 
may require a further demonstration of the durability of the regimen, i.e., a demonstration that the 
drug regimen suppresses the viral load below the level of detectability and keeps it below the level 
of detectability for some period of time. 

As can be appreciated, if a new viral load measurement and tracking technique were 
developed which could reliably detect viral load at levels much lower than the current nadir of 
detectability, the FDA may be able, using the new technique, to much more precisely determine the 
effectiveness of a drug regimen for the purposes of granting approval such that a demonstration of 
the redeeming evalities will no longer be necessary. 

For all of these reasons, it would be highly desirable to provide an improved technique for 
measuring and tracking viral load capable of providing much more precise and reliable estimates of 
the viral load and in particular capable of reducing the nadir of detectability significantly. The 
present invention is directed to this end. 

Summary of the Invention 

In accordance with a first aspect of the invention, a method is provided for determining the 
effectiveness of a therapy, such as an anti-viral therapy, by analyzing biochip output patterns 
generated from biological samples taken at different sampling times from a patient undergoing the 
therapy. In accordance with the method, a viral diffusion curve associated with atherapy of interest 
is generated and each of the output patterns representative of hybridization activity is then mapped 
to coordinates on the viral diffusion curve using fractal filtering. A degree of convergence between 
the mapped coordinates on the viral diffusion curve is determined. Then, a determination is made as 
to whether the therapy of interest has been effective based upon the degree of convergence. 

In an exemplary embodiment, the viral diffusion curve is spatially parameterized such that 
samples map to coordinates near the curve maxima, if the viral load is increasing (i.e., therapy or 
dosage is ineffective). In this manner, any correlation between rate and extent of convergence across 
different patient samples is exploited to provide a quantitative and qualitative estimate of therapy 
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5 effectiveness. 

Also in the exemplary embodiment, the biological sample is a DNA sample. The output 
pattern of the biochip is quantized as a dot spectrogram. The viral diffusion curve is generated by 
inputting parameters representative of viral load studies for the therapy of interest, generating a 
preliminary viral diffusion curve based upon the viral load studies; and then calibrating a degree of 

1 0 directional causality in the preliminary viral diffusion curve to yield the viral diffusion curve. The 
parameters representative of the viral load studies include one or more of baseline viral load (BVL) 
set point measurements at which detection is achieved, BVL at which therapy is recommended and 
viral load markers at which dosage therapy is recommended. The step of generating the 
preliminary viral diffusion curve is performed by selecting a canonical equation representative of 

15 the viral diffusion curve, determining expectation and mean response parameters for use in 
parameterizing the equation selected to represent the viral diffusion curve and parameterizing the 
equation selected to represent the viral diffusion curve to yield the preliminary viral diffusion curve. 

Also, in the exemplary embodiment, each dot spectrogram is mapped to the viral diffusion 
curve using fractal filtering by generating a partitioned iterated fractal system IFS model 

20 representative of the dot spectrogram, determining affine parameters for IFS model, and then 
mapping the dot spectrogram onto the viral diffusion curve using the IFS. Before the dot 
spectrograms is mapped to the viral diffusion curve, the dot spectrograms are interferometrically 
enhanced. After the mapping, any uncertainty in the mapped coordinates is compensated for using 
non-linear information filtering. 

25 In accordance with a second aspect of the invention, a method is provided for determining 

the viral load within a biological sample by analyzing an output pattern of a biochip to which the 
sample is applied. In accordance with the method, a viral diffusion curve associated with a therapy 
of interest is generated and then calibrated using at least two viral load measurements. Then the 
output pattern for the sample is mapped to coordinates on the calibrated viral diffusion curve using 

30 fractal filtering. The viral load is determined from the calibrated viral diffusion curve by 
interpreting the coordinates of the viral diffusion curve. 

Apparatus embodiments are also provided. By exploiting aspects of the invention, disease 
management decisions related to disease progression, therapy and dosage effectiveness may be 
made by tracking the coordinates on the viral diffusion curve as successive DNA-/RNA-based 

35 microarray samples are collected and analyzed. 
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5 Brief Description of the Drawings 

The features, objects, and advantages of the present invention will become more apparent 
from the detailed description set forth below when taken in conjunction with the drawings in which 
like reference characters identify correspondingly throughout and wherein: 

FIG. 1 is a flow chart illustrating an exemplary method for determining the effectiveness of 
10 a viral therapy in accordance with the invention. 

FIG. 2 is a flow chart illustrating an exemplary method for generating Viral Diffusion 
Curves for use with the method of Fig. 1. 

FIG. 3 is a flow chart illustrating an exemplary method for mapping dot spectrograms onto 
Viral Diffusion Curves using fractal filtering for use with the method of Fig. 1. 
15 FIG. 4 is a block diagram illustrating the effect of the fractal filtering of Fig. 3. 

FIG. 5 is a flow chart illustrating an exemplary method for compensating for uncertainty 
using non-linear information filtering for use with the method of Fig. 1. 

Detailed Description of Exemplary Embodiments 

20 With reference to the figures, exemplary method embodiments of invention will now be 

described. The exemplary method will be described primarily with respect to the determination of 
changes in viral loads based upon the output patterns of a hybridized biochip microarray using DNA 
samples, but principles of the invention may be applied to other protein-based samples or to other 
types of output patterns as well. 

25 With reference to FIG. 1, steps will be described for generating viral diffusion curves for 

use is processing DNA biomicroarray output patterns to determine the effectiveness of therapies 
imposed upon a patient providing samples for which the outputs are generated. Then, steps will be 
described for processing the specific output patterns using the VDC's. 

An underlying clinical hypothesis of the exemplary method is that antiviral treatment should 

30 inhibit viral replication and lower an individual's viral load from baseline or suppress rising values. 
A stationary or rising viral load after the introduction of antiviral therapy indicates a lack of 
response to the drug(s) or the development of drug resistance. The VDC exploits the underlying 
hypothesis in part by correlating the rate of disease progression to a sample point value such that a 
change in sample point indicates progression. 

35 The Method 

At step 100, parameters representative of viral load studies for the therapy of interest are 
input. A preliminary viral diffusion curve is generated, at step 102, based upon the viral load 
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studies. The parameters representative of the viral load studies include baseline viral load (BVL) set 
point measurements at which detection is achieved, BVL at which therapy is recommended and 
viral load markers at which dosage therapy is recommended. At step 104, a degree of directional 
causality in the preliminary viral diffusion curve is calibrated to yield the final viral diffusion curve. 

Steps 100 - 104 are performed off-line for setting up the VDC's. These steps need be 
performed only once for a given therapy and for a given set of baseline viral load measurements. 
Thereafter, any number of DNA biomicroarray output patterns may be processed using the VDC's 
to determine the effectiveness of the therapy. Preferably, VDC's are generated for an entire set of 
therapies that may be of interest such that, for any new DNA biomicroarray output pattern, the 
effectiveness of any of the therapies can be quickly determined using the set of VDC's. In general, 
the aforementioned steps need be repeated only to update the VDC's based upon new and different 
baseline viral load studies or if new therapies of interest need to be considered. 

In the following, steps will be summarized for processing DNA biomicroarray output 
patterns using the VDC's to determine whether any therapies of interest represented by the VDC's 
are effective. To determine the effectiveness of therapy at least two samples of DNA to be analyzed 
are collected from a patient, preferably taken some time apart, and biomicroarray patterns are 
generated therefrom. In other cases though, the different samples are collected from different 
patients. 

The output patterns for the DNA biomicroarray are referred to herein as dot spectrograms. 
A dot spectrogram is generated using a DNA biomicroarray for each sample from an N by M DNA 
biomicroarray. An element of the array is an "oxe/'V ofij). An element of the dot spectrogram is a 
hixe[: h(ij). The dot spectrogram is represented by cell amplitudes given by 0(i j) for i: 1 to N, 
and j: 3 to M. 

Dot spectrograms are generated from the samples taken at different times using a 
prefabricated DNA biomicroarray at step 106. The dot spectrograms are interferometrically 
enhanced at step 108. Each dot spectrogram is then mapped to coordinates on the viral diffusion 
curves using fractal filtering at step 110. After the mapping, any uncertainty in the mapped 
coordinates is compensated for at step 112 using non-linear information filtering. 

VDC coordinates are initialized at step 114, then updated in accordance with filtered dot 
spectrograms at step 116. A degree of convergence between the mapped coordinates on the viral 
diffusion curves is then determined at step 118 and a determination is made as to whether the 
therapy of interest has been effective. The determination is based upon whether the degree of 
convergence increases from one DNA sample to another. An increase in degree of convergence is 




WO 00/51055 PCT/USOO/04079 
5 representative of a lack of effectiveness of the therapy of interest. Hence, if the degree of 
convergence decreases, then execution proceeds to step 120, wherein a signal is output indicating 
that the therapy is effective. If the degree of convergence increases, then execution proceeds to step 
122 wherein VDC temporal scale matching is performed. Then a determination is made at step 124 
whether an effectiveness time scale has been exceeded. If exceeded, then a conclusion is drawn that 
1 0 the effectiveness of the viral therapy cannot be established even if more samples are analyzed. If 
not exceeded, then execution returns to step 106 wherein another sample taken from the same 
patient at a latter time is analyzed by repeating steps 106 through 118. 

Viral Load Studies 

15 Viral load studies for therapies of interest are parameterized at step 100 as follows. The 

therapy of interest is selected from a predetermined list of therapies for which viral load studies 
have been performed. Measurements from viral load studies are input for therapy of interest. As 
noted, the viral load measurements include one or more of Baseline Viral Load (BVL) set point 
measurements at which detection is achieved; BVL at which therapy is recommended; and VL 

20 markers at which dosage change recommended. 

Data for the viral load measurements are obtained, for example, from drug qualification 
studies on a minimum include dosages, viral limits as well as time cycles within which an anti-viral 
drug is deemed effective. The data is typically qualified with age/weight outliers and patient history. 
Attribute relevant to this claim is the y, or B VL LOW which corresponds to the lowest detection limit 

25 shown for a therapy to be effective using conventional assays or any other diagnostic means. B VL NP . 
tow denotes the lowest threshold at which viral load is achieved using a nucleotide probe. Using 
interferometric enhancement technique, B VL NP . LOW « B VL^. 

Generation of Viral Diffusion Curves 

30 Referring now to FIG. 2, the viral diffusion curves are generated as follows. An equation is 

selected for representing the VDC at step 200. Expectation (u) and mean response parameters are 
determined at step 202 for use in parameterizing the selected equation. Then the equation selected 
to represent the VDC is parameterized at step 204 to yield a numerical representation of the VDC. 
These steps will now be described in greater detail. 

35 These steps populate a canonical machine representation, denoted as VDC(i, r,y,K) (which 

is a special case of Fokker-Planck equation) to calibrate responses from a viral load detection DNA- 
array based hybridization biomicroarray. 
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i is the index for a diagnostic conditionAherapy of interest, 
T denotes the parameter vector characterizing the VDC, 
X denotes the clinical endpoints vector that indicates detectability 
thresholds for a specific DNA-hybridization array implementation 
10 k correspond to the uncertainty interval estimates. 

An example of an equation selected for representing the VDC is: 



15 |^div(V!P(x)p)+^ , p(*,0) = p°(jt) 



The potential (x) : 9T -> [0, oo) is a smooth function, P > 0 is selected 
constant, and p° (x) is a probability density on <R n . 

20 Preferably, the diffusion potential of the equation and the BVL data are 

such that: 

4>(x) < c . [ BVI^p LOW I BVL LOW ] 
25 the constant c is generally set to 

C = fnumber of amplitude discretization levels] 

[ I og(PCR amplification factor) • avg(oligonucleotides/oxel) ♦ flagging efficiency! ♦ Lbinding efficiency J ) 

30 

Binding efficiency is difficult to quantify analytically for a biomicroarray device 
technology. Hence, for use in the above equation, an estimate of the binding efficiency is 
preferably employed. A binding efficiency of 30% (.3) is appropriate, though other values may 
alternatively used. Depending upon the specific biomicroarray used, the constant c typically ranges 
35 between 0.0001 to 0.5. 
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5 

Expectation and Mean Response Parameters 

The expectation and mean response parameters are then determined at step 202 for use in 
parameterizing the equation selected to represent the VDC. The expectation and mean response 
values are determined by: 1) performing conventional PCR amplification; 2) obtaining calibrated 

10 viral counts from the PCR amplification; 3) determining enhanced and normalized hybridization 
amplitude mean and variance values corresponding to the calibrated viral counts; and 4) matching 
the enhanced and normalized hybridization amplitude mean and variance values. 

Two synthetic amplification techniques (in addition to PCR and any designer tagging) are 
used to achieve VL estimation above the BVL limit set for the exemplary embodiment of the 

15 method, namely (a) readout pre-conditioning, and (b) nonlinear interferometric enhancement. 
Moreover, the expectation match condition implies that: 



Expectation [ Logrt. interferometricallv enhanced image J) 1 vleMtWfifln ^ 0 f interet ti 
Expectation [preconditioned image amplitude] v|expres5ionietofintcrest 

20 

= 1 



Variance matching is done similarly with respect to biomicroarray readout. 
The lower bound of mean response value can be given by: 

25 

Variance \ Log(L inter ferometrical ly enhanced image]) 1 v lmm9km m of iltfcrestl = 1 
Variance [preconditioned image amplitude] v( „ pfes$ion£etofinleT „ 1 

Using the above expression, a conservative lower bound for interferometric enhancement is 
30 estimated for each nucleotide expression of interest. Since the array fabrication device is assumed 
to have (i) an identical oligonucleotide density per oxel and (ii) equal length oligonucleotides, the 
same mean response amplitude can be assumed. If these two assumptions are not met then bounds 
need to be individually calculated and averaged using the above formula. Another assumption is 
that the binding efficiency is statistically independent of the actual oligonucleotide sequence. If this 
35 assumption does not hold for the specific device technology then the binding efficiency should be 
provided as well for each expressed sequence of interest. So the computational analysis method 
uses the analytically derived lower bounds, as computed using the above equation. This is a one- 
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Parameterization of the VDC 

The equation selected to represent the VDC is then parameterized using the expectation and 
mean response values at step 204 to yield a numerical representation of the VDC using 



with constants a, p, y and c « 1. 

This utility of this parameterization is established as follows. The VDC canonical 
representation is based on a variational formulation of the Fokker-Planck. The Fokker-Planck (FP) 
equation, or forward Kolmogorov equation, describes the evolution of the probability density for a 
stochastic process associated with an Ito stochastic differential equation. The exemplary method 
exploits the VDC to model a physical time-dependent phenomena in which randomness plays a 
major role. The specific variant used herein is one for which the drift term is given by the gradient 
of a potential. For a broad class of potentials (that correspond to statistical variability in therapy 
response), a time discrete, iterative variational is constructed whose solutions converge to the 
solution of the Fokker-Planck equation. The time-step is governed by the Wasserstein metric on 
probability measures. In this formulation the dynamics may be regarded as a gradient flux, or a 
steepest descent for the free energy with respect to the Wasserstein metric. This parameterization 
draws from theory of stochastic differential equations: wherein a (normalized) solution to a given 
Fokker-Planck equation represents the probability density for the position (or velocity) of a particle 
whose motion is described by a corresponding lto stochastic differential equation (or Langevin 
equation). The drift coefficient is a gradient. The method exploits "designer conditions" on the drift 
and diffusion coefficients so that the stationary solution of a Fokker-Planck equation satisfies a 
variational principle. It minimizes a certain convex free energy functional over an appropriate 



y(x) = 0o+P,x + P 2 x 2 



subject to constraints 
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admissible class of probability densities. 

A physical analogy is to an optimal control problem which is related to the heating of a 
probe in a kiln. The goal is to control the heating process in such a way that the temperature inside 
the probe follows a certain desired temperature profile. The biomolecular analogy is to seek a 
certain property in the parameterized VDC-namely, an exponential jump in the VDC coordinate 
position for "small linear changes in the viral count". 

This method is in contradistinction to conventional calibration strategies which obtain a 
linear or superlinear shift in quantization parameter for an exponential shift in actual viral count. 

As noted, the form of FP equation chosen is 



where the potential (x) : " -> [0, oo ) is a smooth function, p > 0 is 
selected constant, and p° (x) is a probability density on <R n . The solution p(t,x) is a 
probability density on for almost every fixed time t. So the distribution p(tpc) £ 
0 for almost every (t,x) e (0,qo) X <R D , and 

f r p(/,*)<fr = l 
.for almost every t, e (O^oo). 

It is reasonably assumed that hybridization array device physics for the DNA biomicroarray 
(i.e., corresponding to the potential \y) has an approximately linear response to the nucleotide 
concentration and the response is monotonic with bounded drift. So, 



dp 

dt 




p(x,0) = p°(x) 



where the partition function Z is given by 
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In this model the basis for device physics design is that the potential needs to be modulated 



such that it grows rapidly enough for Z to be finite. This is not achieved by conventional methods. 
However, a technique which does achieve this result is described in co-pending U.S. Patent 
Application xx/xxx,xxx, filed contemporaneously herewith, entitled " Method and Apparatus for 
Analyzing Hybridized DNA Microarray Patterns Using Resonant Interactions Employing Quantum 
Expressor Functions", which is incorporated by reference herein. 

The probability measure p s (x)dx is the unique invariant measure for the Markov random 
field (MRF) fit to the empirical viral load data. 

The method exploits a special dynamical effect to design p . The method restricts the FP 
equation form above to a more restricted case: random walk emulating between critical equilibrium 



To aid in understanding this aspect of the invention, consider the diffusion form 

dt 2 d 2 t 
where D 2 = 7C a 2 
and a is constant. 

A specific VDC shape is parameterized by: 
y(x) = p 0 + p i x + p 2 x 2 
subject to constraints 



with constants a, p, y and e « 1 . 

These constants are set based upon the dynamic range expected for the viral load. Thus, if 



points. 
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5 the viral load is expected to vary only within a factor of 1 0, the constants are set accordingly. If the 
viral load is expected to vary within a greater range, different constants are employed. The actual 
values of the constants also depend upon the particular disease. 
Where the following conditions are met 



The expectation and mean response parameters for use in these equations are derived, as 
described above, from matching the enhanced and normalized hybridization amplitude mean and 
variance that correspond to calibrated viral counts (via classical PCR amplification). 

A distribution represented by the above-equations then satisfies the following form with a 
20 prescribed probability distribution 



10 





2 




15 




Assuming: 




> p > 0, = constant 



25 



and£(t) = £ 0 qy 



such that 



30 
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f(x) = Q.5\/(x)\ 



such that y(-l) < X < y(l). 



The characteristic timescale of response for this system is given by 




arccos 1 - 



B(\/3, 



1,1/3) ayfe~\ 



Ti ~\ 



the successive points must show a motion with characteristic timescale. The VDC is 
designed such that sampling time falls well within characteristic time. 

As noted the actual information used for populating above the parameters is available from 
the following: Baseline viral load (BVL) set point measurements at which detection is achieved; 
BVL at which therapy is recommended; and VL markers at which dosage change recommended. 

The following provides an example of preclinical data that is available to to assist in para- 
meterizing the VDC. 

NOTIONAL VIRAL LOAD MANAGEMENT EXAMPLE 

This is a synthetic example to illustrate how data from clinical studies may be used to 
calibrate the VDC. 

Viral Load analysis studies, using conventional assays, in HIV progression have shown that 
neither gender, age, HCV co-infection, past history of symptomatic HIV-1 infection, duration of 
HIV-1 infection nor risk group are associated with a higher risk of increasing baseline viral load 
(BVL) to the virologic end-point. However, patients with a high (BVL) between 400-6000 copies 
had a 10-fold higher risk of increasing the level of viral load than patients with a BVL below 1500 
copies/ml. Thus, baseline viral load set point measurements provide an important indicator for onset 
of disease. 

Initiation of antiretroviral therapy is generally recommended when the CD4* T-cell count is 
< 600 cells/mL and the viral load level is > 6,000 copies/mL. When the viral load is > 28000 
copies/mL, initiation of therapy is recommended regardless of other laboratory markers and clinical 
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Effective antiretroviral treatment may be measured by changes in plasma HIV RNA levels. 
The ideal end point for effective antiviral therapy is to achieve undetectable levels of virus (< 400 
copies/mL). A decrease in HIV RNA levels of at least 0.5 log suggests effective treatment, while a 
return to pretreatment values (±0.5 log) suggests failure of drug treatment 
10 When HIV RNA levels decline initially but return to pretreatment levels, the loss of therapy 

effectiveness has been associated with the presence of drug-resistant HIV strains. 

The therapy-specific preclinical viral load markers (such as low and high limits in the above 
example) are used to establish actual BVL boundaries for the VDC associated with a particular 
therapy. In this regard, the deterimination of the BVL parameters is disease specific. For example, 
15 in HIV methods such as RT-PCR, bDNA or NASBA are used. Other diseases use other assays. 
Typically, once the parameters of the VDC equation have been set (i.e. constants a, P, y and e), only 
two viral load markers are needed to complete the parameterization of the VDC. This is in contrast 
to previous techniques whereby expensive and laborious techiques are required to determine the 
shape of a viral diffusion curve. The present invention succeeds in using only two viral load 
20 markers in most cases by exploiting the canonical VDC described above which has predefined 
properties and which is predetermined based upon the particular biochip being used. 

Calibration of Viral Diffusion Curves 

25 Referring again to FIG. 1, the directional causality of the VDC is calibrated at step 104 in 

the context of an NIF discussed in greater detail below. At least arbitrarily selected three sample 
points are used execute the NIF calibration computation. The resulting polynomial is used to 
extracting qualitative coherence properties of the system. 

The spectral [0] and temporal coherence [S] is incrementally estimated and computed for 

30 each mutation/oligonucleotide of interest by a NIF forward estimation computation (described 
further below). The two estimates are normalized and convolved to yield cross-correlation function 
over time. The shape index (i.e., curvature) of the minima is used as a measure for directional 
causality. Absence of curvature divergence is used to detect high directional causality in the 
system. 

35 
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5 Sample Collection 

Samples are collected at step 106 subject to a sample point collection separation amount. 
The separation amount for two samples is preferably within half a "drug effectiveness mean time" 
covering 2a population level wherein a denotes the standard deviation in period before which a 
effectiveness for a particular drug is indicated. The following are some general guidelines for 
10 sample preparation for use with the exemplary method: 

- It is important that assay specimen requirements be strictly followed to avoid degradation 
of viral RNA. 

- A baseline should be established for each patient with two specimens drawn two to four 
weeks apart. 

1 5 - Patients should be monitored periodically, every three to four months or more frequently if 

therapy is changed. A viral load level that remains at baseline or a rising level indicates a need for 
change in therapy. 

- Too much significance should not be given to any one viral load result. Only sustained 
increases or decreases of 0.001- 0.01 log [conventional methods typically require a 0.3-0.5 log 

20 change] or more should be considered significant. Biological and technical variation of up to 0.01 
log [typical conventional limit: 0.3 log] is possible. Also, recent immunization, opportunistic 
infections and other conditions may cause transient increases in viral load levels. 

- A new baseline for each patient should be established when changing 
laboratories or methods. 

25 Recommendations for frequency of testing are as follows: 

- establish baseline: 2 measurements, 2 to 4 weeks apart 

- every 3 to 4 months or in conjunction with CD4 + T-cell counts 

- 3 to 4 weeks after initiating or changing antiretroviral therapy 

- shorter intervals as critical decisions are made. 

30 - measurements 2-3 weeks apart to determine a baseline measurement. 

- repeat every 3-6 months thereafter in conjunction with CD4 counts to monitor viral load 
and T-cell count. 

- avoid viral load measurements for 3-4 weeks following an immunization or within one 
month of an infection. 

35 - a new baseline for each patient should be established when changing laboratories or 

methods. 
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The samples are applied to a prefabricated DNA biomicroarray to generate one or more dot 
spectrograms each denoted <X>(i j) for i: 1 to N, and j: 1 to M. The first sample is referred to herein 
as the k=l sample, the second as the k=2 sample, and so on. 



Interferometric Enhancement of the Dot Spectrogram 

10 Each dot spectrogram provided by the DNA biomicroarray is filtered at step 108 to yield 

enhanced dot spectrograms O(k) either by performing a conventional Nucleic Assay Amplification 
or by applying preconditioning and normalization steps as described in the co-pending Patent 
Application having Serial Number 09/253,789, entitled "Method And Apparatus For Analyzing 
Hybridized Biochip Patterns Using Resonance Interactions Employing Quantum Expressor 

15 Functions". The application is incorporated by reference herein, particularly insofar as the 
descriptions of the use of preconditioning and normalization curves are concerned. 



Fractal Filtering 

Each enhanced dot spectrogram is then mapped to the VDC using fractal filtering at step 110 
20 as shown in FIG. 3 by generating a partitioned iterated fractal system 302, determining affine 
parameters for the IFS 304 and then mapping the enhanced dot spectrogram onto the VDC using the 
IFS, step 306. 

The VDC representation models a stochastic process given by 



25 



{°Ay-y D J y*) \yyo) 

(f( x >y) ^ /if 1 (1) , for some 1 <i<m; 



0, 



otherwise; 



for any (x,y)e5R 2 and f <= pOR 2 ) 



30 



An exemplary partitioned iterated fractal system (IFS) model for the system is 
^={0 / =Oi M 7;)}/ = l,2,...,m 

where the affine parameters for the IFS transformation are given by 
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where the D - origin is given by (xj,,^), 
5 the R - origin is given by (x^y^) 

spatial transformation matrix is given by <r f 
the intensity tilting vector is given by r, 
the contrast adjuctment is given by y i% 
the brightness adjustment is given by fi, 

and wherein O represents the enhanced dot spectrogram and wherein \x represents the 
calculated expectation match values 

This IFS model maps the dot spectrogram to a point on the VDC wherein each VDC 
coordinate is denoted by VDC(t,0) such that 
10 W[0,k] -> VDC(k,0) 

Wherein k represents a sample. 

In the above equation, © represent the parameters of the IFS map. 

Thus the output of step 110, is a set of VDC coordinates, identified as VDC(k, 0), with one 
set of coordinates for each enhanced dot spectrogram k = I, 2...n. 

The effect of the steps of FIG. 3 is illustrated in FIG. 4 which shows a set of dot 
spectrograms 450, 451 and 452 and a VDC 454. As illustrated, each dot spectrogram is mapped to 
a point on the VDC. Convergence toward a single point on the VDC implies ineffectiveness of the 
viral therapy. A convergence test is described below. 

Uncertainty Compensation 

With reference to FIG. 5, any uncertainty in the coordinates VDC(k, 0) is compensated 
using Non-linear Information Filtering as follows. Biomicroarray dispersion coefficients, 
hybridization process variability values and empirical variance are determined at step 402. The 
biomicroarray dispersion coefficients, hybridization process variability values and empirical 
variance are then converted at step 404 to parameters for use in NIF. The NIF is then applied at 
step 406 to the VDC coordinates generated at step 106 of FIG. 1. 

Nonlinear information filter (NIF), is a nonlinear variant of the Extended Kai man Filter. A 
30 nonlinear system is considered. Linearizing the state and observation equations, a linear estimator 
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5 which keeps track of total state estimates is provided. The linearized parameters and filter equations 
are expressed in information space. This gives a filter that predicts and estimates information about 
nonlinear state parameters given nonlinear observations and nonlinear system dynamics. 

The information Filter (IF) is essentially a Kalman Filter expressed in terms of measures of 
the amount of about the parameter of interest instead of tracking the states themselves, i.e., track the 
1 0 inverse covariance form of the Kalman filter. Information here is in the Fisher sense, i.e. a measure 
of the amount of information about a parameter present in the observations. 

Uncertainty bars are estimated using NIF algorithm. The parameters depend on 
biomicroarray dispersion coefficients, hybridization process variability and empirical variance 
indicated in the trial studies. 
1 5 One particular advantage of the method of the invention is that it can also be used to capture 

the dispersion from individual to individual, therapy to therapy etc. It is extremely useful and 
enabling to the method in that it can be apriori analytically set to a prechosen value and can be used 
to control the quality of biomicroarray output mapping to VDC coordinates. 

The biomicroarray dispersion coefficients, hybridization process variability values and 
20 empirical variance are determined as follows. Palm generator functions are used to capture 
stochastic variability in hybridization binding efficacy. This method draws upon results in 
stochastic integral geometry and geometric probability theory. 

Geometric measures are constructed to estimate and bound the amplitude wanderings to 
facilitate detection. In particular we seek a measure for each mutation -recognizer centered (MRC-) 
25 hixel that is invariant to local degradation. Measure which can be expressed by multiple integrals of 
the form 

m(Z)=[f(z)dz 

30 where Z denotes the set of mutations of interest. In other words, we determine the function 

f(z) under the condition that m(z) should be invariant with respect to all dispersions Also, up to a 
constant factor, this measure is the only one which is invariant under a group of motions in a plane. 
In principle, we derive deterministic analytical transformations on each MRC-hixel., that map error- 
elliptic dispersion bound defined on <H 2 (the two dimension Euclidean space - i.e., oxel layout) onto 

35 measures defined on 9?. The dispersion bound is given by 
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Such a representation of uniqueness facilitates the rapid decimation of the search space. It is 
implemented using a filter constructed using measure-theoretic arguments. The transformation 
under consideration has its theoretical basis in the Palm Distribution Theory for point processes in 
Euclidean spaces, as well as in a new treatment in the problem of probabilistic description of MRC- 
hixel dispersion generated by a geometrical processes. Latter is reduced to a calculation of 
intensities of point processes. Recall that a point process in some product space E X F is a 
collection of random realizations of that space represented as {(e^ Q, | e, e E, f k e F }. 

The Palm distribution, II of a translation (TJ invariant, finite intensity point process in W is 
defined to the conditional distribution of the process. Its importance is rooted in the fact that it 
provides a complete probabilistic description of a geometrical process. 

In the general form, the Palm distribution can be expressed in terms of a Lebesgue 
factorization of the form 

EpN*= A L N X n 

Where n and A completely and uniquely determine the source distribution P of the 
translation invariant point process. Also, Ep N* denotes the first moment measure of the point 
process and is a probability measure. 

Thus a determination of n and A is needed which can uniquely encode the dispersion and 
amplitude wandering associated with the MRC-hixel. This is achieved by solving a set of equations 
involving Palm Distribution for each hybridization (i.e., mutation of interest). Each hybridization is 
treated as a manifestation of a stochastic point process in 9? 2 

In order to determine n and A we have implemented the following 
measure-theoretic filter: 

Determination of A 

using integral formulae constructed using the marginal density functions 
for the point spread associated with MRC-hixel(i j) 

The oligonucleotide density per oxel p m(iJ) , PGR amplification protocol 
(o ro ), fluorescence binding efficiency (r\ m ) and imaging performance (w J provide 
the continuous probability density function for amplitude wandering in the m-th 



WO 00/51 055 PCT/USOO/04079 

22 

5 MRC-hixe! of interest. Let this distribution be given by (p m(i j V <r m , n m > «0 . 

-The method requires a preset binding dispersion limit to be provided to 
compute A. The second moment to the function 

at SNR = 0 condition is used to provide the bound. 
Determination of II 
1 0 Obtained by solving the inverse problem 

n =&*P 
where 



where t, and x 2 represent the normalized hybridization dispersion limits. 
15 The number are empirically plugged in. The values of 0.1 and 0.7 are appropriate for, 

respectively, signifying loss of 10% - 70% hybridization. Also , 6> denotes the distribution of 
known point process. The form l/(l+exp(p (...))) is employed herein to represent &. 

The biomicroarray dispersion coefficients, hybridization process variability values and 
empirical variance are then converted to parameters at step 304 for use in NIF as follows. 
20 The NIF is represented by: 

Predicted State = f(current state, observation model, information uncertainty, information model) 
Detailed equations are given below. 

In the biomicroarray context, NIF is an information-theoretic filter that predicts and 
estimates information about nonlinear state parameters (quality of observable) given nonlinear 

25 observations (e.g., post hybridization imaging) and nonlinear system dynamics (spatio-temporal 
hybridization degradation). The NIF is expressed in terms of measures of the amount of information 
about the observable (i.e., parameter of interest) instead of tracking the states themselves. It has 
been defined as the inverse covariance of the Kalman filter, where the information is in the Fisher 
sense, i.e, a measure of the amount of information about o/present in the observations Z k where the 

30 Fisher information matrix is the covariance of the score function. 



In a classical sense the biomicroarray output samples can be described by the nonlinear 
discrete-time state transition equation of the form: 



WO 00/51055 PCT/US00/04079 

23 

5 VDC(k) a f{VDC{k-\\<P{k-\\k) + v(k) 

where VDC(k-l) is the state at time instant (k-1), 
<D(k-l) is the input vector (embodied by dosage and/or therapy) 

v(k) is some additive noise; corresponds to the biomicroarray dispersion as computed by the 
Palm Generator functions above. 
1 0 VDC(k) is the state at time k, 

f(k,.,) is the nonlinear state transition function mapping previous state and current input to 
the current state. In this case it is the fractal mapping that provides the VDC coordinate at time k. 

The observations of the state of the system are made according to a non-linear observation 
equation of the form 
1 5 z(k) = h(VDC(k)) + w(k) 

where z(k) is the observation made at time k 

VDC(k) is the state at time k, 

w(k) is some additive observation noise 

and h(.,k) is the current non-linear observation model mapping current state to observations, 
20 i.e., sequence-by-hybridization made at time k, 

v(k) and w(k) are temporally uncorrelated and zero-mean. This is true for the biomicroarray 
in how protocol uncertainties, binding dynamics and hybridization degradation are unrelated and 
additive. The process and observation noises are uncorrelated. 
E[v(i)w T (j)] = 0, Vij. 
25 The dispersion coefficients together define the nonlinear observation model. 

The nonlinear information prediction equation is given by 

y(k\k-\) = y(*|*-l)f (k,VDC(k-\\k-l) 9 u(*-l)) 

Y(k\k-\) = [Vf x (k)Y-> (k - 1 1 k -1) V// (*) + Q(k)]' 1 
The nonlinear estimation equations are given by 
30 y(k\k) = y(k\k-l) + i(k) 

Y(k\k) = Y{k\k-l) + l(k) 
where 

I(k) = W x (k) R- 1 (*)V^(A) 



35 



= Vhl(k) R- 1 (k) [ v{k) + Vh, (k) VDC{k |fc - 1)] 
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5 where ^ 

= | *-!)). 

In this method NIF helps to bind the variability in the VDC coordinate mapping across 
sample to sample so that dosage and therapy effectiveness can be accurately tracked. 

The NIF is then applied to enhanced, fractal-filtered dot spectrogram at step 306 as follows. 
10 States being tracked correspond to post-hybridization dot spectrogram in this method. NIF 
computation as described above specifies the order interval estimate associated with a VDC point. It 
will explain and bound the variability in Viral load estimations for the same patient from laboratory 
to laboratory. 

The NIF also specifies how accurate each VDC coordinate is given the observation model 
1 5 and nucleotide set being analyzed. 



Convergence Testing 

Referring again to FIG. 1, once any uncertainty is compensated, the VDC coordinates are 
renormalized at step 114. The renormalized VDC coordinates are patient specific and therapy 

20 specific. Alternately the coordinates could be virus / nucleotide marker specific. The NIF- 
compensated VDC coordinates are renormalized to the first diagnostic sample point obtained using 
the biomicroarray. Thus a patient can be referenced to any point on the VDC. 

This renormalization step ensures that VDC properties are maintained, notwithstanding 
information uncertainties as indicated by the NIF correction terms. The approach is drawn from 

25 "renormalization-group" approach used for dealing with problems with many scales. In general the 
purpose of renormalization is to eliminate an energy scale, length scale or any other term that could 
produce an effective interaction with arbitrary coupling constants. The strategy is to tackle the 
problem in steps, one step for every length scale. In this method the renormalization methodology is 
abstracted and applied during a posteriori regularization to incorporate information uncertainty and 

30 sample-to-sample variations. 

This is in contradistinction to current viral load measurement calibration methods that either 
generate samples with same protocol and same assumptions of uncertainty or use some constant 
correction term. Both existing approaches skew the viral load readout so that measurements are 
actually accurate only in a limited "information" and "observability" context. This explains the large 

35 variations in readings from different laboratories and technicians for the same patient sample. 



Specifically, we include the dynamic NIF correction function to the 
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5 gradient of the VDC at the sample point normalized in a manner such that when 
the information uncertainty is null, the correction term vanishes. As discussed in 
the above steps, the NIF correction terms is actually derived from the noise 
statistics of the microarray sample. 

<VDC'(k,0)> = VDC(k y 0) + [VMF(y,/)J 
1 0 where VNIF(Y,IX denotes the gradient of nonlinear information prediction function. Under 

perfect observation model this term vanishes. 

Once initialized, the VDC coordinates are then updated at step 116 applying the IFS filter 
W[ ] on k+lth sample, by 

VDC(k+l, 0)<- WfBiomicroarray Output, K+l]; 
1 5 A direction convergence test is next performed at step 1 18 to determine whether the selected 

therapy has been effective. If convergence establishes that the viral load for the patient is moving in 
a direction representative of a lower viral load, then the therapy is deemed effective. The system is 
deemed to be converging toward a lower viral load if and only if: 



20 



KPC(/,)-KDC(/ w ) 
IvDC^-VDCit^l 



VDC^ h -VDC{t k ) 



VDC^ -VDC(f k J 



< 1 fork>2 and j>0 



The above relationships needs to be monotonically persistent for at least two combinations 
of k and j. 

Also, date[k] - date[j] < k * characteristic time, /(in days) 
25 Where k captures the population variability. Typically, k < 1.2 . 

The peak VDC value is determined based on the VDC. The peak amplitude is an artifact of 
the specific parameterization to the Fokker-PIanck equation used in deriving the VDC. It is almost 
always derived independent of the specific sample. 

In connection with step 118, a VDC Shift factor A may be specified at which a dosage 
30 effectiveness decision and/or disease progression decision can be made. The VDC shift factor is 
applied to estimate the VDC curvature traversed between two measurements. 

If the system is deemed to be converging toward a lower viral load, an output signal is 
generated indicating that the therapy of interest is effective at step 120. If not, then the execution 
proceeds to step 122 wherein VDC scale matching is performed. A key assumption underlying this 
35 method is that movement along VDC is significant if and only if the sample points are with a 
constant multiple of temporal scale characterizing the VDC. This does not in any way preclude the 
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5 pharmacological relevance associated with the datapoints. But complete pharmacological 
interpretation of the sample points is outside the scope of this method. The process is assumed to 
be cyclostationary or at a large time scale and two or more sample points have been mapped to 
VDC coordinates. The coordinates are then plugged into an analytic to estimate the empirical cycle 
time (/ ). This is implemented as described in the following sections. 
10 Again the empirical cycle time (? ) is used to establish decision convergence. 



Scale Matching 



15 



20 



25 



Select a forcing function of the form: 
= k. p m cos art 

where k is a constant and m is a small odd integer (m<7). 
The phase space for this dynamical system is represented by: 



. [v^ - ( x V 
x ~ y sin erf m — t= 



k_ 

A+2 



sin oA 



where 



erf m(x) = 



erf(x) 
1 



if x < -N 
if |x|<N 
x>N 



and erf m(.) denotes the error function. 

K is set to 1 ; where o< a, y, co < 1 are refer to constants. 

Let T OTp denote the cycle-time-scale for this empirical system. 

If \og t ( x tm/^) > 1 (instep 10)then we claim that time-scales do not match. 



Time Scale Testing 

Next a determination has been made as to whether an effectiveness timescale has been 
30 exceeded at step 124 by: 

Checking if a time step between successive sampling has exceeded T by determining if 
Time k+I -Timek > T 



such that VDC(k+l, 0)- VDC(k, 0) < q where 
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is set to 0.0001 and wherein 
T is given by 



^. 1 l\ 3(1/3,1/3) aV^l 
10 T arccosh ^ -J 



B(l/3, 1/3) represents the Beta function around the coordinates (1/3,1/3), We can actually use 
all B(l/2i+l,l/2i+l) for i>l and i<7. 

If Time k+ , - Tim^ > T then output signal at step 126 indicating that either 
15 no change in viral load concluded, OR 

therapy deemed ineffective, OR 
dosage deemed suboptimal. 

If Time k+l - Tim^ < T then process another sample by repeating all steps beginning with 
20 Step 4 wherein a dot spectrogram is generated for a new sample. 

If the effectiveness time scale has been exceeded then a signal is output indicating that no 
determination can be as to whether the therapy of interest is effective. If the time scale is not 
exceeded, then execution returns to step 106 for processing another sample. If available, and the 
25 processing steps are repeated. 

Alternative Implementations 

Details regarding a related implementation may be found in co-pending U.S. Patent 
30 Application 09/253,792, entitled "Method and Apparatus for Interpreting Hybridized Bioelectronic 
DNA Microarray Patterns Using Self Scaling Convergent Reverberant Dynamics", and incorporated 
by reference herein. 



35 



The exemplary embodiments have been primarily described with reference to flow charts 
illustrating pertinent features of the embodiments. Each method step also represents a hardware or 
software component for performing the corresponding step. These components are also referred to 
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software component for performing the corresponding step. These components are also referred to 
herein as a "means for" performing the step. It should be appreciated that not all components of a 
complete implementation of a practical system are necessarily illustrated or described in detail. 
Rather, only those components necessary for a thorough understanding of the invention have been 
illustrated and described in detail. Actual implementations may contain more components or, 
depending upon the implementation, may contain fewer components. 

The description of the exemplary embodiments is provided to enable any person skilled in 
the art to make or use the present invention. Various modifications to these embodiments will be 
readily apparent to those skilled in the art and the generic principles defined herein may be applied 
to other embodiments without the use of the inventive faculty. Thus, the invention is not intended 
to be limited to the embodiments shown herein but is to be accorded the widest scope consistent 
with the principles and novel features disclosed herein. 
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CLAIMS 



1 . A method for determining the effectiveness of a therapy by analyzing biochip output 
10 patterns generated from biological samples taken at different sampling times from a patient 
undergoing the therapy, said method comprising the steps of: 

generating a viral diffusion curve associated with a therapy of interest; 
mapping each of the output patterns representative of hybridization activity to respective 
coordinates on the viral diffusion curve using fractal filtering; 
1 5 determining a degree of convergence between the mapped coordinates on the viral diffusion 

curves; and 

determining whether the therapy of interest has been effective based upon the degree of 
convergence from one sample to another. 



20 2. The method of claim 1 wherein the step of generating the viral diffusion curve 

includes the steps of 

inputting parameters representative of viral load studies for the therapy of interest; 
generating a preliminary viral diffusion curve based upon the viral load studies; and 
calibrating a degree of directional causality in the preliminary viral diffusion curve to yield 
25 the viral diffusion curve. 



3 . The method of claim 2 wherein the parameters representative of the viral load studies 
include one or more of baseline viral load (BVL) set point measurements at which detection is 
achieved, BVL at which therapy is recommended and viral load markers at which dosage therapy is 
30 recommended. 



4. The method of claim 2 wherein the step of generating a preliminary viral diffusion 
curve is performed by: 

selecting an equation representative of the viral diffusion curve; 
35 determining expectation and mean response parameters for use in parameterizing the 
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5 equation selected to represent the viral diffusion curve; and 

parameterizing the equation selected to represent the viral diffusion curve to yield the 
preliminary viral diffusion curve. 

5. The method of claim 4 wherein the step of determining expectation and mean 
1 0 response parameters for use in parameterizing the equation selected to represent the viral diffusion 

curve comprises the steps of: 

performing PCR amplification; 

obtaining calibrated viral counts from the PCR amplification; 

determining the enhanced and normalized hybridization amplitude mean and variance values 
1 5 corresponding to the calibrated viral counts; and 

matching the enhanced and normalized hybridization amplitude mean and variance values to 
yield the expectation and mean response parameters. 

6. The method of claim 2 wherein the step of calibrating a degree of directional causality 
20 in the preliminary viral diffusion curve to yield the viral diffusion curve employs nonlinear 

information filtering. 

7. The method of claim 1 wherein the step of mapping each output pattern to the viral 
diffusion curve using fractal filtering comprises the steps of: 

25 generating a partitioned iterated fractal system IFS model representative of the output 

pattern; 

determining affine parameters for IFS model; and 

mapping the output pattern onto the viral diffusion curve using the IFS. 

30 8. The method of claim 1 wherein the step of determining whether the therapy of interest 

has been effective includes the steps of: 

evaluating the degree of convergence against a predetermined threshold of convergence; 
if the degree of convergence does not exceed the predetermined threshold, outputting a signal 
indicating that the therapy of interest has been effective; and 
35 if the degree of convergence exceeds the predetermined threshold, indicating a lack of 
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5 effectiveness of the therapy, then performing viral diffusion curve scale matching on the filtered viral 
diffusion curve and then evaluating whether a predetermined effectiveness time scale has been 
exceeded, and 

if the time scale has been exceeded, outputting a signal indicating the therapy has not been 
effective, and 

10 if the time scale has not been exceeded, then obtaining a new sample and repeating all steps 

using the previously generated viral diffusion curve. 

9. The method of claim 1 further including the step of compensating for uncertainty in 
the filtered viral diffusion curves. 

15 

10. The method of claim 9 wherein the step of compensating for uncertainty in the 
filtered viral diffusion curves is performed using non-linear information filtering. 

1 1 . The method of claim 1 further including the step, performed before the mapping step, 
20 of enhancing the output pattern interferometrically. 

12. The method of claim 11 wherein the step of enhancing the output pattern 
interferometrically comprises the step of performing a nucleic acid amplification. 

25 13. The method of claim 1 wherein the output pattern is a quantized output pattern. 

14. The method of claim 13 wherein the quantized output pattern is a dot spectrogram. 

15. The method of claim 1 wherein the biological sample is selected from a group 
30 consisting of a DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 

(TNA) samples. 



35 



16. A method for determining the viral load within a biological sample by analyzing an 
output pattern of a biochip on which the sample is applied, said method comprising the steps of: 
generating a viral diffusion curve associated with a therapy of interest; 
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calibrating the viral diffusion curve based on at least two viral load measurements; 
mapping the output pattern to coordinates on the calibrated viral diffusion curve using fractal 
filtering; and 

determining the viral load for the sample based upon the viral diffusion curve by interpreting 
the coordinates of the calibrated viral diffusion curve. 

1 7. The method of claim 16 wherein the output pattern is a quantized output pattern. 

1 8. The method of claim 1 7 wherein the quantized output pattern is a dot spectrogram. 

1 9. The method of claim 1 6 wherein the biological sample is a nucleotide-based sample. 

20. The method of claim 1 6 wherein the step of determining the viral load by interpreting 
the coordinates of the VDC is performed by: 

using two or more ground truth measurements. 

21. A system for determining the effectiveness of a therapy by analyzing biochip output 
patterns generated from biological samples taken at different sampling times from a patient 
undergoing the therapy, said system comprising: 

a viral diffusion curve generator for generating a viral diffusion curve associated with a 
therapy of interest; 

a mapping unit for mapping each of the output patterns representative of hybridization 
activity to respective coordinates on the viral diffusion curve using fractal filtering; 

convergence determination unit for determining a degree of convergence between the 
mapped coordinates on the viral diffusion curves for determining whether the therapy of interest has 
been effective based upon the degree of convergence from one sample to another. 

22. The system of claim 21 wherein the viral diffusion curve generator includes 

an input unit for inputting parameters representative of viral load studies for the therapy of 
interest; 

a preliminary viral diffusion curve generator for generating a preliminary viral diffusion 
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5 curve based upon the viral load studies; and 

a calibration unit for calibrating a degree of directional causality in the preliminary viral 
diffusion curve to yield the viral diffusion curve. 

23 . The system of claim 2 1 wherein the parameters representative of the viral load studies 
10 include one or more of baseline viral load (BVL) set point measurements at which detection is 
achieved, BVL at which therapy is recommended and viral load markers at which dosage therapy is 
recommended. 



15 operates by 

selecting an equation representative of the viral diffusion curve; 

determining expectation and mean response parameters for use in parameterizing the 
equation selected to represent the viral diffusion curve; and 

parameterizing the equation selected to represent the viral diffusion curve to yield the 
20 preliminary viral diffusion curve. 

25. The system of claim 24 wherein the preliminary viral diffusion curve generator 
operates to determine expectation and mean response parameters for use in parameterizing the 
equation selected to represent the viral diffusion curve by: 

25 performing PCR amplification; 

obtaining calibrated viral counts from the PCR amplification; 

determining the enhanced and normalized hybridization amplitude mean and variance values 
corresponding to the calibrated viral counts; and 

matching the enhanced and normalized hybridization amplitude mean and variance values to 
30 yield the expectation and mean response parameters. 

26. The system of claim 22 wherein the calibration unit employs nonlinear information 
filtering. 



24. The method of claim 22 wherein the preliminary viral diffusion curve generator 



35 



27. 



The system of claim 21 wherein the mapping unit operates by: 
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5 



generating a partitioned iterated fractal system IFS model representative of the output 



pattern; 



determining affine parameters for IFS model; and 

mapping the output pattern onto the viral diffusion curve using the IFS. 



10 



28. The system of claim 2 1 wherein the convergence determination unit operates by: 



evaluating the degree of convergence against a predetermined threshold of convergence; 

if the degree of convergence does not exceed the predetermined threshold, outputting a signal 
indicating that the therapy of interest has been effective; and 

if the degree of convergence exceeds the predetermined threshold, for indicating a lack of 
15 effectiveness of the therapy and for performing viral diffusion curve scale matching on the filtered 
viral diffusion curve and then evaluating whether a predetermined effectiveness time scale has been 
exceeded, and 

if the time scale has been exceeded, outputting a signal indicating the therapy has not been 
effective, and 

20 if the time scale has not been exceeded, obtaining a new sample and repeating all steps using 

the previously generated viral diffusion curve. 

29. The system of claim 2 1 further including a compensation unit for compensating for 
uncertainty in the filtered viral diffusion curves. 



30. The system of claim 29 wherein the compensation unit employs non-linear 
information filtering. 

3 1 . The system of claim 3 1 further including an interferometric enhancement unit for 
30 enhancing the output pattern interferometrically. 

32. The system of claim 3 1 wherein the interferometric enhancement unit operates to 
perform a nucleic acid amplification. 



25 



35 



33. The system of claim 3 1 wherein the output pattern is a quantized output pattern. 
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34. The system of claim 33 wherein the quantized output pattern is a dot spectrogram. 

35. The system of claim 31 wherein the biological sample is selected from a group 
consisting of a DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 

10 (TN A) samples. 

36. A system for determining the viral load within a biological sample by analyzing an 
output pattern of a biochip on which the sample is applied, said system comprising: 

a viral diffusion curve generator for generating a viral diffusion curve associated with a 
1 5 therapy of interest; 

a calibration unit for calibrating the viral diffusion curve based on at least two viral load 
measurements; 

a mapping unit for mapping the output pattern to coordinates on the calibrated viral diffusion 
curve using fractal filtering; and 
20 a viral load determination unit for determining the viral load for the sample based upon the 

viral diffusion curve by interpreting the coordinates of the calibrated viral diffusion curve. 

37. The system of claim 36 wherein the output pattern is a quantized output pattern. 
25 38. The system of claim 37 wherein the quantized output pattern is a dot spectrogram. 

39. The system of claim 36 wherein the biological sample is selected from a group 
consisting of a DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 
(TNA) samples. 

30 

40. The system of claim 36 wherein the determination unit operates using measured 
ground truth states. 
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